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Abstract 

Results are presented for the leading two-loop contributions of 0(cttct s ) to the masses 
and mixing effects in the Higgs sector of the MSSM with complex parameters. They are 
obtained in the Feynman-diagrammatic approach using on-shell renormalization. The 
full dependence on all complex phases is taken into account. The renormalization of 
the appropriate contributions of the Higgs-boson sector and the scalar top and bottom 
sector is discussed. Our numerical analysis for the lightest MSSM Higgs-boson mass is 
based on the new two- loop corrections, supplemented by the full one- loop result. The 
corrections induced by the phase variation in the scalar top sector are enhanced by 
the two-loop contributions. We find that the corresponding shift in M/ il can amount 
to 5 GeV. 
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1 Introduction 



The Higgs sector of the Minimal Supersymmetric Standard Model (MSSM) with two scalar 
doublets accommodates five physical Higgs bosons. In lowest order these are the light and 
heavy CP-even h and H, the CP-odd A, and the charged Higgs bosons H . Higher-order con- 
tributions yield large corrections to the masses and couplings, and also induce CP-violation 
leading to mixing between h, H and A in the case of general complex SUSY breaking pa- 
rameters. 

For the MSSM with real parameters (rMSSM) the status of higher-order corrections to 
the masses and mixing angles in the Higgs sector is quite advanced [1-10]. In the case of the 
MSSM with complex parameters (cMSSM), the first more general investigations [11] were 
followed by evaluations in the effective potential approach [12] and with the renormalizat ion- 
group-improved one-loop effective potential method [13, 14]. These results have been re- 
stricted to the corrections arising from the (s)fermion sector and some leading logarithmic 
corrections from the gaugino sectoiQ. Within the Feynman diagrammatic (FD) approach the 
one-loop leading m\ corrections have been evaluated in Ref. [15]. Most recently a full one- 
loop calculation in the FD approach was presented [16] (further discussions on the effect of 
complex phases on Higgs boson masses can be found in Ref. [17]) and implemented in the pro- 
gram FeynHiggs [2,7,16,18], which is publicly available. Another public code, CPsuperH [19], 
is based on the renormalization-group-improved effective potential approach [13,14]. 

In this letter we improve our diagrammatic one-loop calculation [16] by providing the 
leading 0(a t a s ) corrections of the Higgs-boson masses and mixings in the cMSSM obtained 
in the FD approach. Technically we calculate and renormalize the Higgs-boson self energies 
taking into account the general complex parameters of the appropriate part of the colored 
sector of the cMSSM. We provide numerical examples for the lightest cMSSM Higgs-boson 
mass and discuss the dependence on the phases in the scalar top sector and on the gluino mass 
parameter. The results presented in this paper will be included in the code FeynHiggs [20]. 



2 The Higgs-boson sector of the cMSSM 

With the two Higgs doublets of the cMSSM decomposed in the following way, 
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1 The two-loop results of [10] can in principle also be taken over to the cMSSM. However, no explicit 
evaluation or computer code based on these results exists. 
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where the ellipses stand for higher powers in the Higgs-boson fields. In Eq. (J2J) the tadpoles 
appear as the coefficients of the linear terms, and the bilinear terms contain the neutral and 
charged mass matrices M^ xx and M^i^i. Tadpoles and mass matrices are conveniently 
rewritten in terms of the physical components h, H, A, H and the Goldstone components 
G, G ± . Details about the tadpole coefficients and the mass matrices can be found in Ref. [16]. 

Eq. flTJ introduces a possible new phase £ between the two Higgs doublets. The potential 
Vh contains the real soft breaking terms m\ and m\ (with m\ 



m\ 



m 2 = m\ + |/i 



and the generally complex soft breaking parameter m\ 2 , entering the mass matrices and 
tadpoles in Eq. (T5]). With the help of a Peccei-Quinn transformation [21], \i and m\ 2 can 
be redefined [22] such that the complex phase of m\ 2 vanishes. In the following we will 



therefore treat real parameter, i.e. 
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Together with the requirement 



that the minimum of Vh is located at v 1 and v 2 , all tadpoles are zero at lowest order. 

Investigating the Higgs potential beyond the tree level, renormalization has to be applied 
to the mass matrices and the tadpoles, introducing counterterms according to the loop 
expansion up to second order, 
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where the mass matrices and tadpoles are obtained from those in Eq. (J2J) by a rotation 
to the physical states. The leading 0(a t a s ) contributions to the Higgs-boson self-energies 
are obtained in the limit of vanishing gauge couplings and neglecting the dependence on 
the external momentum. The bottom Yukawa coupling, appearing in the charged Higgs- 
boson self-energy, is also neglected. Since for the leading two-loop terms the Goldstone- 
boson parts of Eqs. ([3]) and (j4j) do not contribute (see also the discussion in Ref. [16]), mass 
renormalization at the two-loop level reduces to the counterterms 
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The renormalized Higgs-boson self-energies, denoted as Ejj(p 2 ) with i,j = h, H, A, H ± , 
are expanded into a one- loop and a two- loop part, 



S w (p a ) = EgV) + s!?(0). 
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The complete one-loop part has been obtained in Ref. [16], and the two-loop part is eval- 
uated at vanishing external momentum, as explained above. The renormalized two-loop 
self-energies 
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(21 

involve the unrenormalized Higgs-boson self-energies Ej - (0), containing the one-loop sub- 
renormalization, and the counterterms of Eq. ([6]). 

The entries Sm^ = h,H,A) of the counterterm matrix in Eq. are not all 

2(2) (2) 

independent, but can be expressed in terms of SM^J and ST! . As explained e.g. in Ref. [2], 
M^ ± and Tj are the only independent parameters in the Higgs potential that have to be 
renormalized for the evaluation of the 0(a>ta> s ) terms. Correspondingly, it is sufficient to 
impose renormalization conditions for the tadpoles and for the charged Higgs-boson mass: 

• The tadpoles are fixed by the requirement that the minimum of the Higgs potential is 
not shifted, yielding at the two-loop level 
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The mass square of the charged Higgs boson, M^±, is fixed by an on-shell condition 
yielding the following counterterm at the two-loop level: 
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With these, the counterterms for the neutral mass matrix are now determined in the following 

way, 
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We have used s x = sin(x), c x = cos(x) as abbreviations. The angle a diagonalizes the 
mass matrix at tree-level, and denote the H and h tadpoles, respectively, in the limit 
of a — > 0, which are the tadpoles in Eq. (j2J). 

The calculation of the unrenormalized self-energies and tadpoles at 0(at<y s ) requires the 
evaluation of genuine two-loop diagrams and one-loop graphs with counterterm insertions. 
Example diagrams for the neutral Higgs-boson self-energies are depicted in Fig. [TJ and for 
the charged Higgs boson in Fig. [2j Examples for the tadpole diagrams are displayed in Fig. [31 
The complete set of contributing Feynman diagrams has been generated with the program 
FeynArts [23]; tensor reduction and the evaluation of traces was done with support of the 
programs OneCalc and TwoCalc [24], yielding analytic expressions in terms of the scalar 
one-loop functions A , B [25] and two-loop vacuum integrals [26]. The numerical evaluation 
was performed with the help of the program LoopTools [27]. 
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Figure 1: Examples of generic two- loop diagrams and diagrams with counterterm insertion 
for the Higgs-boson self-energies (0 = h, H, A; k, I — 1, 2). 
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Figure 2: Examples of generic two- loop diagrams and diagrams with counterterm insertion 
for the charged Higgs-boson self-energy k = 1, 2). 



Figure 3: Examples of generic two- loop diagrams and diagrams with counterterm insertion 
for the Higgs-boson tadpoles, (0 = h, H, A; i,j, k = 1, 2). 



The renormalized self-energies determine the dressed propagators of the Higgs fields, 
from which masses and mixing properties at higher order are derived. The self-energies have 
an impact on the location of the poles and thus on the Higgs particle masses, which are in 
general different from their tree-level values. Only the charged Higgs boson mass M^± is 
not shifted, owing to the on-shell renormalization condition (TlTj) . 

The non-diagonal self-energies are responsible for mixing in the neutral Higgs system. In 
the presence of complex parameters all three neutral CP eigenstates h, H, A can mix. The 
3x3 propagator matrix, A^^(p 2 ), is obtained by inverting the renormalized irreducible 
two-point function, 




'j 
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where 



(19) 



(m\-t hh {f) -t hH (P 2 ) -t hA {P 2 ) \ 

M n (p 2 )=\ -£ hH (p 2 ) m 2 H -t HH {p 2 ) -± HA (p 2 ) . (20) 
V -± hA {p 2 ) -t HA {p 2 ) m\-± AA {p 2 )) 

The masses of the three Higgs-boson mass eigenstates, h\, h 2 , h 3 , ordered according to 
Mh t < Mh 2 < Mh 3 , are given by the real parts of the poles of AhH A (p 2 ) or, equivalently, 
of the roots of the determinant of the two-point vertex function, det[ThH A (p 2 )} = 0. The 
quantities m h ,m H ,m A in Eq. (|20!) are the masses of h,H,A at the tree-level, respectively. 



3 The colored sector of the cMSSM 

For the evaluation of the 0{ata s ) two-loop contributions to the tadpoles and self-energies, a 
renormalization of the one-loop contributions from the scalar top (t) and bottom (b) sector 
is needed, giving rise to the counterterms for one- loop subrenormalization (see Figs. [EH2D- 
The bilinear part of the t and b Lagrangian, 

% mass = " HI 4) Mf (£) - (bl P R ) M~ b , (21) 

contains the stop and sbottom mass matrices and M^, given by 

/ M? + ml + Mlc 20 (T* - Q q s 2 ) m Q X* 

q V ™ q X q M\ R + m\ + M 2 c 2p Q q s 2 1 

with 

X q = A q — ji*K, k = {cot (5, tan/5} for q = t,b. (23) 

Q q and T g 3 denote charge and isospin of q, and A q is the trilinear soft-breaking parameter. 
The mass matrix can be diagonalized with the help of a unitary transformation XJ q , which 
can be parametrized by a mixing angle 9 q and a phase ip q , 

t _ (m% \ TT fUti„ ^7ff 10 ^ / cos 6^ e^sinfl; 



T") - — T T - lVT - T JT — { 91 " TT _ ^<?ii W 912 _ a ^"<z 

9 " 9 9 9 " V m 2 J ' U? " \U q21 U,J ~ sin 6 q cos ^ 

(24) 

The mass eigenvalues depend only on \X q \. 

Taking into account complex phases, the renormalization in the t sector is somewhat 
more involved than in the case of real parameters [2,3,28]. In the cMSSM the t sector is 
described in terms of five real parameters (where we assume that \x and tan j3 are defined via 
other sectors): the real soft SUSY-breaking parameters M\ and M 2 ^, the absolute value and 
complex phase of the trilinear coupling, A t = \A t \e %ipA t, and the top Yukawa coupling At that 
can be chosen to be real. Instead of the quantities M 2 , M 2 ^ and X t , in the on-shell scheme 
applied in this paper we choose the on-shell squark masses m? , m~ 2 and the top-quark mass 
m t as independent parameters. 

The following renormalization conditions are imposed: 
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(i) The top-quark mass is defined on-shell, yielding the one- loop counterterm 5m t : 

Sm t = -m t (Re£f (m 2 ) + ReSf {m 2 t ) + 2ReSf {m 2 t )) , (25) 

referring to the Lorentz decomposition of the self energy E t 

=^_S f L (A; 2 )+ ^+Sf (/c 2 ) + m f Ef (A; 2 ) + m tl ^ s (k 2 ) (26) 

into a left-handed, a right-handed, a scalar and a pseudoscalar part, S^, Ef , and 
£f 5 , respectively. Re denotes the real part with respect to contributions from the loop 
integral, but leaves the complex couplings unaffected. 

(ii) The stop masses are also determined via on-shell conditions [2,28], yielding 

5m\ = ReE t ~..(m|) with i = 1, 2 . (27) 

'iii) The third condition affects the stop mixing angle and phase, or equivalently, the A t 
parameter. Rewriting the squark mass matrix in terms of the mass eigenvalues and 
the mixing angle and phase using Eq. ( l24l ). 

/ cos 2 Ozm 2 - +sin 2 #fm~ e l<p tsm9i cosOAmZ —m 2 ~)\ 
M- = ( 1 2 o 1 2 ) (28) 

* \e~ m sin 9 1 cos 9 i(m 2 — m 2 - ) sin 9pn~ + cos 2 9pn~ I ' v 7 

yields the counterterm matrix 5M t - by introducing counterterms Sm 2 ^ , <5w 2 2 for the 
masses and <5# f ~, for the angles. For the case of the real MSSM one obtains the 
counterterm for the mixing angle, 

rMSSM: (m| - m|) S9 { = [U f< JM t -Uj-] 12 = SY U (29) 

for which the following renormalization condition has been used [3,28]: 

rMSSM: 5Y~ t = ~ [ReE t ~ 12 (m| ) + ReE t ~ 12 (m|)] . (30) 

Generalizing Eq. (|29|) to the complex case, we obtain 

SYi = [U^M,-Ut] 12 = {ml-m\)e i ^{59 i + ism9 iC os9 i 5ip~) , (31) 
and impose, as a generalization of Eq. (|30|) . the condition 

5Y t = \ [ReS fi2 (m|) + ReE fl2 (m|)] , (32) 

which now corresponds to two separate conditions for the real and imaginary part, or 
for 59i and Sip^, respectively. 

We adopt a scheme where \A t \ and tpA t are chosen as independent parameters. The two 
sets of parameters # f ~, ipi and \A t \, ipA t are mutually related via Eq. f[2"2"j) and Eq. f[2"Bl . 
The off-diagonal entries of the corresponding counterterm matrices yield 

(A* t -» cot (3) Smt + m t 5A* t = U? u U il2 {5m\ - 8m\ ) + U~ t22 5Y t + U il2 5Y t * . (33) 
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As a result, we obtain for 5\A t \ and 5<fA t 



5\A t \ = —Re[e^K t ] , (34) 
m t 

8<pm = — l rnM^ At K t ] , (35) 

m t \At\ 

with 

K t = -(A* t -ficotfiSrrn + C/^OK - 5m|) + t^t^,^ + U^SY? . 

In the scalar bottom sector, we also encounter five real parameters (with \x and tan/3 
defined via other sectors): the real soft-breaking parameters M\ and M~ , the absolute 
value and phase of the trilinear coupling A&, and the bottom Yukawa coupling A& that can 
be chosen to be real (for the set of corrections presented in this paper A& does not enter, 
as explained above). SU(2) invariance requires the "left-handed" soft-breaking parameters 
in the stop and the sbottom sector to be identical (denoted as M|). In the evaluation of 
the 0(a t a s ) contributions to the Higgs-boson self-energies, the counterterms of the sbottom 
sector appear only in the self-energy of the charged Higgs boson. In our approximation, 
where the 6-quark mass is neglected, bi and bn do not mix, and b^ decouples and does not 
contribute. The charged Higgs-boson self-energy thus depends only on a single parameter 
of the sbottom sector, which can be chosen as the squark mass vn- hL . The parameter m~ bL 

should be regarded simply as the upper left entry in the b mass matrix, not as a physical 
b pole mass (see also Ref. [3,28]). By means of SU(2) invariance, the corresponding mass 
counterterm is already determined: 

5m\ L = \U- tll \Hm\ + \U k f5ml - U^JY - U il2 W t JY; - 2m t 5m t . (36) 

With the set of renormalization constants determined in Eqs. (1231) . (I2~?j) . (|34|) . (135]) and 
(136]) the counterterms for the diagonal and non-diagonal (s) quark self-energies as well as for 
all Higgs-boson-(s) quark vertices are at our disposal for the one-loop subrenormalization. An 
explicit list of the counterterms will be provided in a detailed forthcoming publication [20]. 

Finally, at 0(a t a s ) gluinos appear as virtual particles at the two-loop level (hence, no 
renormalization is needed). The corresponding soft-breaking gluino mass parameter M3 is 
in general complex, 

M 3 = \M 3 \e i<Ps (with the gluino mass m~ g = |M 3 |) . (37) 

The phase can be absorbed by a redefinition of the gluino Majorana spinor such that it 
appears only in the gluino couplings, but not in the mass term. 

4 Numerical results 

We illustrate the effects of the two-loop contributions in terms of the mass of the lightest 
neutral Higgs boson, M^, evaluated on the basis of Eq. (120]) with the entries from Eq. (]7j). 
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The results for physical observables are affected only by certain combinations of the complex 
phases. In particular, the 0(a t a s ) corrections presented in this paper depend only on the 
combinations [22,29] 

fiA t (m^)* and A t M 3 * . (38) 

As discussed above, we have transformed our parameters such that the complex phase of 
m\ 2 vanishes. Therefore our two- loop results depend on the phases of the parameters A t , fi 
and M 3 , which we denote as <pA t , anc ^ Pgi respectively. We do not consider the variation 
of complex phases that enter only via one-loop contributions. 

In the context of a detailed phenomenological analysis of the cMSSM parameter space the 
existing constraints on CP-violating parameters from experimental bounds for the electric 
dipole moments (EDMs) [30, 31] are of interest. While SM contributions enter only at 
the three-loop level, due to its complex phases the cMSSM can contribute to the EDMs 
already at one-loop order. The complex phases appearing in the cMSSM are experimentally 
constrained by their contribution to the EDMs of heavy quarks [32], of the electron and 
the neutron (see [33,34] and references therein), and of deuterium [35]. One finds that 
in particular the phase ip^ is tightly constrained (in the convention where the phase of the 
gaugino mass parameter M2 is set to zero). The bounds on the phases of the third-generation 
trilinear couplings, on the other hand, are much weaker. 

Since the complex phases appear in our two-loop result only in the combinations given 
in Eq. fl38|) . we can conveniently choose p^ = 0, so that in our numerical analysis only tpA t 
and <fg are varied. In order to illustrate the possible effects of complex phases we will show 
below results for p g , <pA t varied over the full parameter range. 

Our numerical analysis has been performed for the following set of parameters (if not 
indicated differently): 

Msusy = 1000 GeV, \A t \ = \A b \ = \A T \ = 1000 GeV, <p Ab = <p Ar = 0, 

fj. = 1000 GeV, M 2 = 500 GeV, M x = (5s 2 J / (3c 2 w ) M 2 , m g = 1000 GeV, 

M H ± = 500 GeV, tan/3 = 10, m t = 174.3 GeV . (39) 

-^susy denotes the diagonal soft SUSY-breaking parameters in the sfermion mass matrices 
that are chosen to be equal to each other. We do not consider higher values of tan /3, which 
in general enhance the SUSY contributions to the EDMs. 

We first discuss the dependence of M hl on the phase in the scalar top sector. Since 
the leading one-loop result in the limit M H ± 3> Mz depends only on the absolute value 
\X t \ = \A t — fi* I tan/5| (implying that only the combination <pA t + f/i enters, in accordance 
with Eq. fl38|) ). it is useful to analyze the dependence of the result on ip At as well as on ipx t @ 
In Fig. H]we show the lightest Higgs-boson mass as a function of <p At (left) and of ip x t (right) 
for \X t \ = 1.5 TeV (upper row) and \X t \ = 2.5 TeV (lower row). \A t \ is chosen such that for 
vanishing phases it is equal in the left and right plot of each row. A variation of <pA t for fixed 
li and tan (3 changes the absolute value of X t and thus the masses of the scalar top quarks. 
Changing <px t , on the other hand, leaves the masses of the scalar tops invariant (see Sect. [3]), 
but changes A t . Therefore, in the right plots the t masses are constant (m^ = 770 GeV and 
mi 2 = 1210 GeV). We compare in Fig. H]the one-loop result for (dotted line) with the 
new result that includes the 0(a t a s ) contributions (solid line). 

2 It should be noted that the variation of Mh x with ipx t can be substantial for small values of M H ± [17]. 
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Figure 4: In the left column is given as a function of ipA t for \At\ = 1.6 TeV (upper 
plot) and \A t \ = 2.6 TeV (lower plot). In the right column is given as a function of ipx t 
for \X t \ = 1.5 TeV (upper plot) and \X t \ = 2.5 TeV (lower plot). The other parameters are 
as given in Eq. (I59"j) and ipg = 0. The one-loop results (dashed line) are compared with the 
results including the 0(a t a s ) corrections (solid line). 



Fig. |4] shows that the two-loop contributions lead to a reduction of M/ ll of ~ 20 GeV, in 
accordance with the results in the rMSSM [1,2]. The dependence on the complex phases ipA t 
and (fx t is much more pronounced in the two-loop result than in the one-loop case. This 
fact can easily be understood from the discussion above: for the relatively large M H ± chosen 
in Fig. H] the one-loop result is dominated by contributions involving only the absolute value 
\X t \. Therefore the dependence of the one-loop result on tpx t is very weak (right column 
of Fig. H]), while the dependence on ip^ t arises to good approximation only from its effect 
on \X t \. At the two-loop level, on the other hand, the contributions with internal gluinos 
depend on the phase of (A t M£), see Eq. (1381) . This induces an asymmetry of the leading 
corrections to M/ ll with respect to X t (see Refs. [2,4] for a discussion in the rMSSM). 

The impact of the phases (fA t , <Px t is obviously enhanced for larger values of \A t \ and 
\X t \. In Fig. [5] we show the dependence of on \A t \ (left) and \X t \ (right) for ipA t ,fx t — 
0,7r/2,7r, respectively. Concerning the dependence on \X t \, at the one-loop level the results 
are indistinguishable for the three values of *fix t i in agreement with the one- loop results in 
Fig. HI Varying tpA t , on the other hand, results in a shift of the position of the maximum 
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Figure 5: as a function of \A t \ (left plot) and \X t \ (right plot) for (pxt^At — 0,7r/2,7r. 

The other parameters are as given in Eq. (I3"§j) and (fg = 0. The one- loop results are compared 
with the results including the 0(a t a s ) corrections. 
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Figure 6: M hl at (9(a t a s ) as a function of (/^ for M SU sy = 500 GeV (left) and M SUS y 
1000 GeV (right) with M H ± = 150, 500 GeV, tan/3 = 3, 5, 10 and tp At = 0. 



of M/ ll in the one-loop result (left plot). At the two-loop level, the position and size of the 
maximum value of M/^ is significantly affected both by (fA t an d <fix t , i n accordance with the 
discussion above. 

We now investigate the dependence of on the phase of the gluino mass parameter, 
keeping the phase of A t fixed at (f At = 0. This means that only the second term in Eq. (1381) 
is affected by the phase variation, while so far we had studied the combined effect of both 
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Figure 7: Mh t at 0(a t a s ) as a function of m g for Msusy = 500 GeV (left) and Msusy = 
1000 GeV (right) with ip g = 0, tt/2, tt and <p At = 0. 

terms in Eq. fl38l) . Fig. [6] displays the variation of M hl with ip g for M SU sy = 500 GeV (left) 
and Msusy = 1000 GeV (right). M H ± is set to 150,500 GeV, and tan/5 = 3,5,10. The 
dependence on the gluino phase (for tp^ t — 0) is relatively weak for the set of parameters 
chosen in Fig. [61 yielding shifts in M/ ll below ~ 2 GeV. For larger M#± the dependence is 
slightly stronger than for small M H ± values. In all cases a minimum of Mh 1 is reached for 
(f g = tt. Larger effects of the phase of the gluino mass parameter than the ones shown in 
the example of Fig. [6] would occur for larger values of \A t \, as a consequence of Eq. (|38|) . 

In Fig. [7Jthe result for is shown as a function of m g for tp g = 0, tt/2, it (and <fA t = 0). 
Msusy is set to 500 GeV (left) and 1000 GeV (right). The phase dependence is strongest 
around the thresholds m g = m$ — m t and m g = m 4 ~ 2 — m t . For the chosen set of parameters 
the thresholds correspond to rrig = 177 GeV (not shown) and rrig = 487 GeV for Msusy = 
500 GeV, and to m g = 760 GeV and m g = 915 GeV for M SUS y = 1000 GeV. The change in 
M/ ll induced by the phase variation can amount up to 4 GeV in the threshold area for the 
parameters chosen in Fig. [71 

5 Conclusions 

We have presented results for the leading 0(a t a s ) contributions to the dressed Higgs-boson 
propagators in the MSSM with complex parameters, obtained in the Feynman-diagrammatic 
approach using an on-shell type renormalization scheme. In the Higgs sector a two-loop 
renormalization has to be carried out for the mass of the charged Higgs boson and the three 
tadpoles. The renormalization of the scalar top and bottom sector at the one-loop level 
involves a renormalization of the complex phase (fA t ■ 

Concerning the explicit numerical results we have focused on the lightest Higgs-boson 
mass, Mh x . This is of interest in view of the current exclusion bounds [36,37] and possible 
high-precision measurements of the properties of a light Higgs boson at the next generation of 
colliders [38-40]. The 0(ata s ) corrections yield a large downward shift in M^, in accordance 
with the well-known result from the rMSSM. We find that the impact of the complex phases 
ip A t and <fx t is significantly enhanced by the two-loop contributions, which is a consequence 
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in particular of diagrams involving internal gluinos. We find that varying the complex phases 
of the scalar top sector and of the gluino mass parameter can induce shifts in of up to 
~ 5 GeV even in cases where the one-loop result shows hardly any dependence on the phases. 
The result for M hl for <^ t , ipg ^ 0, ir is found to lie in the intervals given by ±|A t |, ±|M 3 |. 
The effects of the complex phases of the 0(a t a s ) corrections can also be enhanced in the 
threshold region where the gluino mass is approximately equal to the sum of the top-quark 
mass and the mass of one of the scalar top quarks. 

The new results of 0(atCt s ) will be implemented into the Fortran code FeynHiggs [2,7, 
16,18]. A detailed description, including a comparison of different renormalization schemes 
for the scalar top sector and a more elaborate discussion of Higgs-boson masses and mixings 
will be presented in a forthcoming publication [20] , as well as a comparison with the results 
based on the renormalization-group improved effective-potential approach [14,19]. 
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